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Abstract 

Habitat degradation through anthropogenic development is a key driver of biodi- 
versity loss. One way to compensate losses is "biodiversity offsetting" (wherein 
biodiversity impacted is "replaced" through restoration elsewhere). A challenge in 
implementing offsets, which has received scant attention in the literature, is the 
accurate determination of residual biodiversity losses. We explore this challenge 
for offsetting gas extraction in the Ustyurt Plateau, Uzbekistan. Our goal was to 
determine the landscape extent of habitat impacts, particularly how the footprint 
of "linear" infrastructure (i.e. roads, pipelines), often disregarded in compensa- 
tion calculations, compares with "hub" infrastructure (i.e. extraction facilities). 
We measured vegetation cover and plant species richness using the line-intercept 
method, along transects running from infrastructure/control sites outward for 
500 m, accounting for wind direction to identify dust deposition impacts. Find- 
ings from 24 transects were extrapolated to the broader plateau by mapping total 
landscape infrastructure network using GPS data and satellite imagery. Vegetation 
cover and species richness were significantly lower at development sites than con- 
trols. These differences disappeared within 25 m of the edge of the area physically 
occupied by infrastructure. The current habitat footprint of gas infrastructure is 
220 ± 19 km 2 across the Ustyurt (total ~ 100,000 km 2 ), 37 ± 6% of which is lin- 
ear infrastructure. Vegetation impacts diminish rapidly with increasing distance 
from infrastructure, and localized dust deposition does not conspicuously extend 
the disturbance footprint. Habitat losses from gas extraction infrastructure cover 
0.2% of the study area, but this reflects directly eliminated vegetation only. 
Impacts upon fauna pose a more difficult determination, as these require 
accounting for behavioral and demographic responses to disturbance by elusive 
mammals, including threatened species. This study demonstrates that impacts of 
linear infrastructure in regions such as the Ustyurt should be accounted for not 
just with respect to development sites but also associated transportation and 
delivery routes. 



Introduction 

Land-use change is a key driver of biodiversity decline 
through impacts on habitat availability (Mace et al. 2005). 
The extractive hydrocarbon, metal and mineral industries 
are among the most locally damaging forms of anthropo- 
genic disturbance in many areas (Baillie et al. 2004). The 
negative impacts of extractive activities could be at least 
partially compensated for through "biodiversity offsets," 
which endeavor to ensure "no net loss" of biodiversity 
alongside development by compensating for unavoidable 



ecological losses (ten Kate et al. 2004; Bull et al. 2013a). 
However, formal and comprehensive quantification of the 
landscape impacts of infrastructure in the relevant environ- 
ment will be essential if the potential of offsetting is to be 
realized (Quintero and Mathur 2011). Accurate quantifica- 
tion of the residual impacts from development, which 
biodiversity offsets aim to compensate for, is often over- 
looked in the offset literature. Semiarid habitats may be par- 
ticularly vulnerable to disturbance from industrial activities 
(Lovich and Bainbridge 1999), recovering poorly if at all 
(Fiori and Martin 2003) and with full suites of constituent 
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species persisting only in remnant undisturbed patches 
(Rapport and Whitford 1999). Given that dryland (i.e. arid 
and semi-arid) biomes cover approximately 41% of the 
world's land surface and support over 38% of the global 
human population (Millennium Ecosystem Assessment 
2005), it is crucial that responses to industrial disturbance 
in these habitats are better understood and managed. 

In semiarid environments, habitats heavily disturbed by 
anthropogenic activities often have lower species richness 
compared with undisturbed areas (Simmers and Gala- 
towitsch 2010). Even low-intensity and small-scale distur- 
bances can have immediate and persistent effects (Forbes 
et al. 2001), with slow recovery times (Cui et al. 2009). The 
impacts of infrastructure upon vegetation in a variety of 
habitats - and of roads in particular, which are ubiquitous 
in regions of industrial extractive activity - have been well 
studied (Forman and Alexander 1998; Trombulak and Fris- 
sell 2000; Coffin 2007). Irrespective of the sector that uses 
them, roads can impact ecosystems in numerous ways: via 
non-native species brought in by vehicles (Gelbard and Bel- 
nap 2003); nitrous oxide and other pollutants produced by 
vehicles (Gadsdon and Power 2009); and by subdividing 
populations and forming physical barriers to dispersal that 
alter demographics (Forman and Alexander 1998). Thus, 
species richness (Lee et al. 2012) and total vegetation cover 
(Fiori and Martin 2003) may decrease in response to road 
construction. The magnitude of these effects will vary with 
distance from the actual disturbance. For example, in Cali- 
fornia, Gelbard and Harrison (2003) found that plant cover 
was significantly lower within 10 m of roads and that 
native species richness was impacted even 1 km away. Con- 
versely, roads may actually enhance overall species richness 
(Zeng et al. 2011), for example, if associated disturbances 
provide favorable microsites for vegetation establishment 
(Boeken and Shachak 1994; Brown and Schoknecht 2001). 

Natural resource extraction has had well-documented 
impacts on wildlife as well as vegetation (E&P Forum/ 
UNEP 1997; Epstein and Selber 2002; Kumpula et al. 
2011; OGP/IPIECA 2011). Habitat disturbance from 
industrial infrastructure can affect wildlife both spatially 
and temporally, for instance altering breeding patterns of 
birds (Walker et al. 2007) and grazing patterns of herbi- 
vores, leading to increased usage and pressure on sur- 
rounding undisturbed habitats (Vistnes and Nellemann 
2007). Noise associated with transport and other activities 
along infrastructure may disrupt the use of habitat by 
faunal species (Rabanal et al. 2010), and the physical 
obstruction caused by pipelines can alter animal move- 
ment across landscapes (Curatolo and Murphy 1986; Dyer 
et al. 2002). In this manner, the spatial impacts of distur- 
bance on the mean species abundance of vertebrate 
mammals in a range of habitats have been assessed to 
extend up to 5 km from infrastructure (Bem'tez-Lopez 
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et al. 2010). Vegetation responses may likewise be impor- 
tant to wildlife, altering resource availability and habitat 
structure. Further indirect effects may come from 
increased human use of the area that follows develop- 
ment, including water extraction, natural resource use, 
and hunting (Thibault and Blaney 2003). Finally, as with 
roads, industrial activity can be associated with the spread 
of invasive alien species (E&P Forum/UNEP 1997), 
although this has not been a conspicuous problem for the 
case study region (J. Bull, pers. obs.). 

Many of the ecological impacts from industry can, in 
principle, be mitigated through biodiversity offsetting 
("offsets"), but offsets are beset with implementation 
challenges, including the difficulty of fully quantifying 
biodiversity losses (Bull et al. 2013a). Quantification of 
the scale of disturbance impacts for offset projects often 
focuses on the impacts of development sites or "hubs," 
whereas the disturbance from the construction and opera- 
tion of linear infrastructure such as roads that link "hubs" 
is only recently being treated as something that could also 
be compensated for through offset mechanisms (Quintero 
and Mathur 2011). Further, if linear infrastructure has a 
comparable impact by area to hub infrastructure across a 
landscape, which to the authors' knowledge has not 
before been explored for offsets, then the scale of offsets 
required to achieve "no net loss" of biodiversity could be 
greater than currently thought. 

Biodiversity and natural gas extraction in 
semiarid Uzbekistan 

Central Asian countries are predominately semi-arid and 
constitute an important proportion of the global semiarid 
biome, 16% of which is found in Asia. In Uzbekistan, 
over 99% of the country by area is arid and semi-arid 
(White and Nackoney 2003). Uzbekistan contains an esti- 
mated 27,000 species (USAID 2001; UNDP 2010). The 
transboundary Ustyurt Plateau, which Uzbekistan shares 
with Kazakhstan, covers ~ 100,000 km 2 of Uzbekistan and 
contains 271 recorded vascular plant species (Gintzburger 
et al. 2011), several of which are on the IUCN Red List 
(Esipov and Shomurodov 2011). The Ustyurt is also 
home to vertebrate species of high conservation priority, 
most notably the critically endangered saiga antelope 
Saiga tatarka tatarka (Mallon 2008). The saiga is the 
flagship conservation target for the Ustyurt, and its distri- 
bution and abundance are driven by both vegetation 
cover (Singh et al. 2010a) and sensitivity to anthropo- 
genic presence (Singh et al. 2010b). 

The Ustyurt is not only valuable from an ecological 
perspective, but also for the subterranean resources it 
contains. Uzbekistan is among the top 20 global exporters 
of natural gas (Effimoff 2000; EIA 2012), and this industry 
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is expanding (Dorian 2006). Karakalpakstan, the adminis- 
trative region containing the Uzbekistan part of the 
Ustyurt, is one of the country's main growth areas for gas 
extraction, with more than 3bn USD of foreign investment 
currently agreed for exploration and extraction activities 
(UNDP 2010; EIA 2012). Development of the Ustyurt for 
hydrocarbon production since the Soviet era has resulted 
in extensive infrastructure growth, including exploration 
and extraction sites, pipelines, and a substantial network 
of roads (UNDP 2010). To date, as far as we are aware, 
there has been no quantitative evaluation of the impact of 
roads, or other components of natural gas infrastructure 
(pipelines, extractive facilities), upon biodiversity on the 
Ustyurt. 

Quantifying disturbance on the Ustyurt 
Plateau 

There are currently plans to design an offset policy in the 
Uzbekistan Ustyurt, to compensate for the ecological losses 
from development and conserve the regional saiga popula- 
tion, as part of a broader biodiversity conservation initia- 
tive for the Uzbek oil and gas sector (UNDP 2010; Bull 
et al. 2013b). Offset projects on the Ustyurt could be 
improved by the provision of sound scientific knowledge of 
the impacts of natural gas activity at the landscape level, 
including whether linear infrastructure impacts are impor- 
tant enough to necessitate inclusion in offset calculations. 
The types of impact we explore here - direct impacts on 
vegetation from industrial infrastructure, in an area of gen- 
eral conservation concern but without legal protection - 
could be considered to be unavoidable, residual impacts of 
the natural gas industry. These direct impacts are therefore 
appropriate for biodiversity offsetting (Quintero and 
Mathur 201 1), which is an idea at the forefront of consider- 
ations for conservation in the Ustyurt (Bull et al. 2013b). 
Condition and area-based assessments of vegetation losses 
and gains form the basis of a number of existing offset poli- 
cies (Quetier and Lavorel 2011), hence our focus on mea- 
suring impacts on vegetation quality (condition) and cover 
(area). We aim not only to quantify the local impacts of 
infrastructure on Ustyurt habitat, but also the landscape- 
scale habitat footprint of the oil and gas sector in a way that 
is relevant for the development of an appropriate biodiver- 
sity offset policy. 

On the Ustyurt, the vast majority of roads created and 
used by the oil and gas industry are unpaved and formed 
from repeated vehicle travel rather than formal clearing 
and construction. Most roads are therefore temporary or 
seasonal, although arterial routes between established 
extraction facilities are large and effectively permanent. In 
addition to damage from direct impaction and clearing 
by vehicles, dust deposition could further impact plateau 
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vegetation, as has been suggested by observations of dust 
cloud movement and settling following vehicle passage 
(Gintzburger et al. 2011). Dust can affect key physiologi- 
cal processes such as photosynthesis, respiration, and 
transpiration (Farmer 1993). In arid environments, dust 
can coat leaf surfaces, altering their radiation balance 
(Grantz et al. 2003), and can alter nutrient cycling 
through effects on soil bacteria and fungi, potentially 
costly in a nutrient-limited environment (Forbes et al. 
2001). Dust in this region may be particularly harmful to 
organisms due to the very high content of pesticide resi- 
dues and heavy metals that have resulted from drying of 
the Aral Sea (Micklin 2007). Strong winds are characteris- 
tic of this continental landmass, and so the dominant 
wind direction and subsequent dust deposition from dis- 
turbance sources have been suggested as potential drivers 
of vegetation response to disturbance on the Ustyurt 
(Micklin 2007). 

In this study, we quantify habitat impacts of oil and gas 
extraction activities by first investigating the spatial extent 
of localized disturbance from infrastructure, measured in 
both plant species richness and vegetation cover. These two 
metrics together provide a broad measure of habitat condi- 
tion in this region (Opp 2005), have successfully been used 
in previous studies investigating the impacts of disturbance 
(e.g. Fiori and Martin 2003; Gelbard and Harrison 2003; 
Simmers and Galatowitsch 2010; Lee et al. 2012), and pro- 
vide a baseline for understanding the state of plant commu- 
nities at a range of scales. We conducted field surveys 
throughout the region to determine whether: ( 1 ) there is a 
reduction in vegetation species richness and ground cover 
in disturbed relative to undisturbed (control) sites; (2) spe- 
cies richness and vegetation cover increase with distance 
from disturbance and reach background levels within 
500 m from the disturbance; and (3) effects of disturbance 
persist over greater distances in areas downwind of the dis- 
turbance, according to dominant wind direction. 

We use the results from field surveys to estimate the 
spatial extent of the local disturbance footprint of natural 
gas infrastructure across the entire plateau. By plotting 
the plateau-scale network of gas infrastructure, we extrap- 
olate from localized measurements on the ground to the 
magnitude of the existing natural gas disturbance foot- 
print in its entirety. This provides a measure of impact 
for vegetation that can inform ongoing development of a 
biodiversity offset policy for the region. 

Materials and Methods 

Vegetation impact transects 

Following a pilot expedition in September and October 
2011 to test sampling protocols and arrange access, an 
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18-day field expedition to the Ustyurt Plateau, Uzbeki- Disturbances were designated as either linear (road or 

stan, in May and June 2012, led by the United Nations pipeline) or point ("hubs" e.g. gas extraction sites). In 

Development Programme (UNDP), was undertaken for each case, transects ran from the center of the disturbance 

formal data collection. Transects employing the line-inter- outward for 500 m. These were perpendicular to the dis- 

cept method (Canfield 1941) were used to gather data on turbance for linear sources and at a randomly directed 

plant species richness (total number of species) and vege- angle for hubs. The size and type of each primary distur- 

tation cover. Surveys were carried out at randomly bance were recorded, and then transects were arrayed in a 

selected locations within eight regions across the Ustyurt "spine and rib" formation, with increased sampling effort 

Plateau. The eight survey areas were selected by local close to the disturbance source to detect any fine-scale 

experts to reflect the heterogeneous nature of Ustyurt bio- vegetation responses (Fig. 2; Angold 1997; Lee et al. 

diversity (Fig. 1; Esipov and Shomurodov 2011). Within 2012). The line-intercept measurement method was 

each designated survey area, we completed as many tran- applied over the 20 m length of each rib. To apply the 

sects as possible depending on time available, ensuring line-intercept method, we laid down a 20 m length of 

that they were spaced > 500 m apart. Transects originated rope marked with 5-cm intervals. For each plant overlaid 

at the center of components of oil and gas infrastructure by the rope, we recorded the length of rope intercepted 

and included extraction sites, pipelines, and unpaved by the plant when viewed from directly above. The num- 

roads known to be primarily used by oil and gas compa- ber of different species of plant was also noted, meaning 

nies. For the purposes of this study, these major sources that for each rib, we could calculate both total vegetation 

of disturbance were deemed "primary disturbances." percentage intercept (a proxy for percent cover) and spe- 

Transects were also undertaken at control sites without cies richness - species richness and percentage cover were 

any of the above three types of infrastructure present for the response variables used in analyses, 
a minimum of 2 km. In total, 24 transects were com- There is a network of small (< 3 m width) unpaved 

pleted: 17 disturbed sites and seven control sites. tracks across the Ustyurt not primarily affiliated with the 




Figure 1. Geography of Central Asia, Uzbekistan, and the Ustyurt plateau. Circles on the detailed map are survey locations, with transects 
roughly equally distributed across the eight locations. Maps created using Garmin "BaseCamp"™ software. 
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Figure 2. Design of spine and rib transects: the main 500-m spine 
transect originates from the center of disturbance with 20-m rib 
transects bisecting it at set intervals: 25-m intervals between 0 and 
100 m (where 0 m is the center of disturbance), 50-m intervals 
between 100 and 300 m, and 100-m intervals between 300 and 
BOO m, giving increased sampling effort closer to the disturbance. The 
line-intercept method was used to collect species richness and 
vegetation cover data along each rib. Secondary disturbances (i.e. 
small tracks) were also recorded. 



extractive industry that spanned control and disturbed 
transects; these were labeled as secondary disturbances 
and noted and measured in surveys. Binary data on the 
presence/absence and width of secondary disturbances 
within 10 m and 50 m of ribs were recorded. 

Due to the limited size of the dataset, we could not test 
the effects of each disturbance type - roads, pipelines, 
and extraction sites - separately, but we note that patterns 
of vegetation impact and recovery appeared visually simi- 
lar for each disturbance type. 



Methodological considerations and 
compromises 

It was important to ensure that sites were representative 
of the Ustyurt landscape as a whole, to allow scalable 
inferences about the impacts of infrastructure. Habitat on 
the Ustyurt is generally classified by local botanists into 
different key species associations, for example the wide- 
spread Anabasis sp. - Artemisia sp. - Salsola sp. associa- 
tion. Nevertheless, different associations are comparable 
in terms of both percentage cover and species diversity 
(Gintzburger et al. 2003; Esipov and Shomurodov 2011). 
Surveys were thus completed throughout the Uzbek 
Ustyurt and within a range of association types, but dif- 
ferent association types were not separated in the analyses 
reported here. 

Survey areas were chosen such that, in general, each 
transect remained within one association type. In order to 
minimize observer site selection bias as far as possible for 
linear disturbances, the start point was selected by walk- 
ing in a fixed direction along the road/pipeline for a 
pre-agreed amount of time from the point at which the 
survey vehicle stopped (5 min), or 500 m from the end 
of the previous transect if more than one was being com- 
pleted in a sequence. The point reached after this time 
was the transect start point. 

Statistical analysis 

All statistical analyses were conducting using linear 
mixed-effects models to allow inclusion of both fixed and 
random effects (Bolker et al. 2009) in the "R" statistical 
package (R Development Core Team 2011). All models 
had species richness or vegetation cover as response vari- 
ables. As species richness was in the form of count data, a 
Poisson error structure was specified for analyses. For 
percentage cover, data were arcsine-square root-trans- 
formed so as to follow a Gaussian rather than a binomial 
distribution in analyses. Maximal models included the 
variables site type (control vs. disturbed), distance from 
disturbance, disturbance width (either <3 m or >10 m), 
presence and width of secondary disturbances, and tran- 
sect direction (as a proxy for wind direction). Wind 
direction data were obtained for 2009/10 from the hydro- 
meteorological stations in Jaslyk (central Ustyurt) and in 
Muynak (eastern Ustyurt). As data from both sites clearly 
showed dominant easterly and northerly wind direction, 
and local wind data were not available for the transect 
sites visited, we assumed that this was the dominant wind 
direction throughout. 

The transect design resulted in "ribs" being pseudore- 
plicated and nested within "transect." To account for this, 
"transect" was fitted as a random effect within linear 
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mixed-effects models. Maximal models including all 
potential explanatory variables were simplified through 
stepwise deletion of highest order nonsignificant terms, 
and model comparison using ANOVA (Crawley 2007). 
Models used in the analyses were as follows: 

lmer(species richness ~ control or disturbed sites x dis- 
tance from disturbance + disturbance width <3 or 
>10 m + (1\ transect), family = poisson); 
and, after arcsine-square root transformation of percent- 
age cover data: 

Imer (percentage cover ~ control or disturbed x distance 
from disturbance + (l\transect). 

P-values for general linear mixed-effects models involv- 
ing percentage cover as the response variable were pro- 
duced using Markov chain Monte Carlo (MCMC) 
sampling with 10,000 iterations, following Bolker et al. 
(2009). The package "languageR" and function "pvals.fnc" 
(Baayen 2011) were used with the lme4 package (Bates 
et al. 2012) to run MCMC and obtain P-values for analy- 
ses involving percentage cover data. 

Landscape footprint of oil and gas activities 

Detailed official records for the spatial extent of oil and gas 
infrastructure on the Ustyurt were not available. Therefore, 
we used two methods to map the spatial distribution of 
key oil and gas infrastructure. First, during two field expe- 
ditions to the Ustyurt (the 2012 main expedition and the 
2011 pilot expedition), a GPS track log was kept of all oil 
and gas infrastructure (tracks, pipelines, roads, railways, 
compressor stations, wells) encountered throughout the 
Plateau. These data were used to create an infrastructure 
map, extrapolated as necessary to join up known pipeline 
routes, and with the addition of the coordinates of addi- 
tional known gas extraction facilities (UNDP 2010). 

Secondly, satellite imagery was used to map all linear 
infrastructure on the Ustyurt Plateau clearly visible at 
70 km altitude (Google Earth 2012). At this altitude, infra- 
structure known to be present from field observations was 
clearly visible, validating the method; it was assumed that 
other infrastructure similarly visible in the images, but not 
visited during the fieldwork, was also attributable to natu- 
ral gas activity given the lack of other major infrastructural 
activities in the area. Lower altitude images, although 
more detailed, would have included less distinct infra- 
structure relating to general use of the area by herders or 
local traffic. It is noted that this method was not also used 
to map hub infrastructure: it was found difficult to distin- 
guish with any confidence between hub gas extraction sites 
and other unrelated development hubs, such as settle- 
ments, takyrs (salt pans), and even water wells. 

Hub disturbances consisted of the six major known gas 
extraction facilities: two which were visited during the 



fieldwork (Shakpakty and Aqsholaq) and four unnamed 
additional sites mapped by the UNDP (2010) but inacces- 
sible for visitation. Limited information was available on 
the spatial footprint of gas extraction facilities, and those 
that we visited in the field could not be directly surveyed. 
However, a recent environmental impact assessment com- 
pleted for the proposed Surgil gas development in Uzbeki- 
stan provides detailed area calculations for a major gas 
facility on the Ustyurt (Mott-Macdonald 2012). These cal- 
culations were taken as reflective of the spatial footprint 
of a standard major gas facility in the region and were 
used as representative data to calculate the footprints of 
the six hub sites. The major gas compressor plants at Ja- 
slyk, Karakalpakia ,and Kubla-na-Ustyurt were excluded, 
because although they represent major infrastructure, they 
are based within substantial settlements that would have 
existed in any case. Therefore, their footprint could not 
be attributed directly to industrial gas activities. 

The physical footprint (in km 2 ) occupied by each infra- 
structure type was calculated separately for both the field- 
based and satellite-based mapping techniques. In both 
cases, we calculated the footprint using mapped lengths 
with our field-measured widths of infrastructure compo- 
nents to obtain the overall area. The outcomes of these two 
different mapping approaches were compared qualitatively. 
Averaging the footprints obtained from both approaches 
produced an estimated total footprint of linear infrastruc- 
ture, and the difference was treated as a simplistic measure 
of uncertainty. This was then summed with the total hub 
infrastructure footprint, to give an estimated total foot- 
print of oil- and gas-related infrastructure on the Ustyurt. 

Results 

Vegetation impacts of infrastructure 

Oil- and gas-associated disturbance had a negative effect 
on species richness (z = —6.2, P < 0.001) and vegetation 
cover (t = —4.7, P < 0.001) compared with control sites 
(Table 1 and Fig. 3). 

Outside of actual infrastructure, species richness and 
cover increased to baseline levels exhibited by control sites 
(Fig. 4). Species richness and cover were only significantly 
different to baseline levels at the site of disturbance itself 
(0 m); at 25-500 m from disturbance, species richness and 
cover were not significantly different from baseline levels. 
Our sampling design did not anticipate complete attenua- 
tion of disturbance-associated vegetation differences within 
25 m of infrastructure and was therefore not set up specifi- 
cally to test for patterns between 0 and 25 m distance. 

There were no significant directional impacts upon 
either richness or cover, as would be expected if wind 
deposition of dust or exhaust pollutants was an important 
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Figure 3. Boxplots showing differences in vegetation responses at 
the 0-m "rib" transect of control and disturbed sites. Differences in 
both species richness (z = -6.2, P< 0.001) and cover (f = -4.7, 
P < 0.001) are significant. 



si/1 




0 25 50 75 100 200 300 500 
Distance (m) 

Figure 4. Interaction plots for richness and cover with distance. 
Hollow points represent disturbed sites, solid points represent 
controls. Graphs produced using "Sciplot" with 95% confidence 
intervals displayed. 

factor for vegetation. Species richness was higher for small 
(< 3 m width) and large (> 10 m width) linear 
disturbances than it was for medium (3-10 m) linear dis- 
turbances (z = 2.4, P < 0.05 and z = 2.8, P < 0.01, respec- 
tively). 

Footprint of infrastructure at landscape 
level 

The two mapping approaches gave visually comparable 
estimates of linear infrastructure presence (Fig. 5). In the 




Figure 5. Mapping of spatial extent of gas infrastructure. (A) 
mapped using GPS data, where lines = roads, tracks, and pipelines, 
circles = settlements with gas infrastructure, triangles = known 
facilities used or created by the gas industry. Black lines = mapped 
using GPS data, gray lines = known infrastructure not mapped using 
GPS data; (B) linear infrastructure mapped using Google Earth (2012). 
White lines represent linear infrastructure - roads, railways, pipelines 
- associated with the oil and gas industry. Note that the large white 
region near the center of this figure corresponds to a dried-up saline 
lake. 



NE of the plateau, infrastructure that was missed using 
Google Earth was recorded using GPS on the ground. 
This casts some doubt on the ability of the remotely 
sensed approach to capture all key components of infra- 
structure in an investigation of this type. Conversely, 
however, in the SE of the plateau, the Google Earth map 
shows linear infrastructure that was not surveyed during 
the expedition. Although this is simply because it was not 
possible to cover the plateau in its entirety during the 
expedition, the potential value of using the Google Earth 
approach to remotely map infrastructure is clear. While 
neither approach is comprehensive, the two arguably 
complement each other. The field-collected data suggested 
a current linear infrastructure footprint =100 km 2 ; the 
Google Earth method suggested a footprint = 63 km 2 . 

It was not possible to survey hub extraction facilities 
for health and safety reasons, but we could see from 
visual inspection that vegetation was cleared in a similarly 
comprehensive manner within hub infrastructure as it 
was for linear infrastructure. With an assumption of 
immediate attenuation of disturbance effects, as conserva- 
tively documented for linear infrastructure, the footprint 
of extraction sites on vegetation can be calculated as the 
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Table 1. Overview of vegetation response to primary disturbances. The distance at which vegetation is significantly different to baseline levels 
(those at the 500-m sampling point) is shown. Number of transect sites = 17, number of control sites = 7. 



Effect of primary disturbance 



Species richness 
Response 



Percentage cover 
Response 



Distance at which vegetation is significantly different to baseline levels 



Species richness 



Percentage cover 



Distance (m) 



Response 



Distance (m) 



Response 



z= -6.17 

P = 6.87e-10*** 



t= -4.68 

P = 0.0000*** 



z = -4.66 
P= 3.14e-6*** 



f = -5.2 

P = 0.0000*** 



***Statistically significant 



actual ground surface area occupied by extraction facili- 
ties. Extrapolating from the 23 km 2 proposed Surgil 
facility (based on data provided in Mott-Macdonald 2012; 
Table 2), yields an estimate of a 138 km 2 footprint. 

Summing the estimates for linear and hub infrastructure, 
and using the difference in linear infrastructure estimates 
from the two different methods as a measure of uncer- 
tainty, suggests a total current footprint of 220 ±19 km 2 . 
Again, this is in the context of a plateau that covers 
~ 100,000 km 2 of Uzbekistan. Under this rough estimate, 
linear infrastructure constitutes 37 ± 6% of the total oil 
and gas infrastructure footprint. 

Discussion 

Local impact of disturbance 

This study shows that natural gas infrastructure, both 
linear and hub, has a substantial effect on local species 
richness and vegetation cover in a semiarid ecosystem. 
But, contrary to expectation, effects are effectively limited 
to the footprint of physical infrastructure itself, at least 
in so far as our sampling scheme was equipped to mea- 
sure. Future sampling that focuses on the zone between 
0 and 25 meters distant from oil and gas infrastructure 
is key to deriving a more precise estimate. Our results 
did not support the idea that dust deposition signifi- 
cantly impacts the local cover or richness of vegetation. 
Future work could further quantify dust deposition 
including its more subtle or widespread effects (Goossens 
and Rajot 2008) and investigate whether directional 
effects on vegetation metrics vary spatially across the Us- 
tyurt. In conjunction, these two findings suggest that the 
condition of vegetation was only reduced within the 
direct footprint of infrastructure, and there it was effec- 
tively reduced to zero: making the development of a con- 
dition-area type biodiversity offset metric, in this case, 
straightforward. 

This study does not assess all impacts associated with 
oil and gas development. For instance, the industry cer- 
tainly abstracts a large amount of water and may result in 
the keeping of additional livestock on the plateau, both of 



which may impact upon regional biodiversity; equally, 
other elements of infrastructure are thought to impact 
specific components of biodiversity, such as power lines 
causing bird mortality (Mott-Macdonald 2012). However, 
it does quantify the limited effects of infrastructure on 
vegetation as one component of a potential biodiversity 
offset policy in the region. Our work adds to current 
knowledge of disturbance effects from infrastructure on 
plant communities in climatically severe areas. These 
include Arctic tundra (Kemper and Macdonald 2009), 
steppe environments (Fiori and Martin 2003), and other 
semiarid desert regions (Simmers and Galatowitsch 2010). 
Oil and gas exploration is burgeoning in these habitat 
types; consequently, providing accurate assessments of the 
spatial scale of infrastructure disturbance is essential if 
national governments intend to ensure that any negative 
ecological impacts associated with the oil and gas industry 
are compensated for. But, as previous work strongly sug- 
gests (Bemtez-Lopez et al. 2010), the spatial scale of local 
disturbance is likely to extend beyond the direct physical 
footprint for infrastructure, when considering other com- 
ponents of the ecosystem such as faunal species. Hence, 
the area calculated here is a minimum estimate. Taking 
other effects of infrastructure into account is also likely to 
increase the relative importance of linear infrastructure in 
calculating overall compensation requirements, because of 
the known disruptive effects of anthropogenic activity on 
wildlife (Bemtez-Lopez et al. 2010), particularly migratory 
species such as the saiga antelope (e.g. Singh et al. 
2010b). Further work could include investigation of the 
extent to which disturbance effects can be seen in other 
taxonomic groups on the Ustyurt such as invertebrates, 
reptiles, mammals, and birds, and whether the spatial 
scales of responses differ. 

There were often observed to be relatively undisturbed 
"humps" between tire tracks, and there is the possibility 
that water pooling in tire ruts may create favorable estab- 
lishment sites in arid ecosystems (Boeken and Shachak 
1994; Briones et al. 1998; Brooks and Lair 2005). Further- 
more, in some cases and especially for smaller, less busy 
roads, it is possible that roads had not been used for 
some time, allowing some minor colonization by vegeta- 
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Table 2. Breakdown of estimated contributions to total footprint for 
proposed Surgil Gas Extraction Facility (based on data in Mott-Mac- 
donald 2012). Assumed dimensions and calculations by the authors 
are given in italics in the "Area" column. Assumed road widths, 
extraction radii, pipeline widths, and railway widths were based on 
filed observations. Values in regular font come from the report itself. 
Connecting ("linear") infrastructure on site is included in the site 
("hub") footprint. 



Component 


Subcomponent 


Area (m 2 ) 


Main facilities 






Ustyurt Gas 


Main 


970,100 


Chemical 


Storage 


27,000 


Complex 


Wastewater 


240,000 




pond 




Settlement 


Area 


850,000 




Road to 


5 km by 10 m = 50,000 




settlement 






On-site camp 


? (treated as zero) 


Gas extraction 


Wells 


133 wells (radius 






of 170 m)= 12,075,340 




Gas gathering 


6 stations (radius 




stations 


of 170 m) = 544,752 


Complex Gas 


Expansion of 


? (treated as zero) 


Treatment Unit 


existing area 




Water 


Water treatment 


? (treated as zero) 




plant 




Connecting 






infrastructure 






Roads 


Access (UGCC) 


9 ha = 90,000 




Gas gathering 


6 roads, 3 km by 




stations 


25 m = 450,000 




Gas wells 


133 wells, average distance 






A. V m m^rt \A/i/~ifh c\f 

H- Mil, / UCJU VVILILII \JI 






10 m = 5,320,000 


Pipelines 


Gas 


115 km by 13 m = 1,495,000 




Gas sales 


9 km by 13 m = 117,000 




Gas booster 


? (treated as zero) 




station 






Water 


27 km by 13 m = 351,000 


Railway 


Access 


7 km by 30 m = 210,000 


Electricity 


Transmission 


18 km by 10 m = 180,000 




line (road) 




Total 








Main facilities 


14,757,192 




Connecting 


8,213,000 




infrastructure 






Sum 


22,970,192 



tion. This may explain why cover and richness were not 
always zero within the footprint of disturbances. 

Scaling up to landscape level 

Our use of mapping tools in conjunction with local-scale 
quantification of direct infrastructure impacts allowed esti- 
mation of the disturbance footprint from infrastructure for 
the whole region. This estimate is based on a number of 



broad assumptions and is therefore an order-of-magnitude 
estimate only. The two approaches we used were comple- 
mentary: mapping using GPS data ground-truths 
infrastructure sites, while mapping using Google Earth 
captures all infrastructure of a certain size. Furthermore, 
Google Earth images present a single snapshot of the 
Ustyurt (without clear dates for images in the version used: 
Google Earth 2012), and so does not reveal the seasonality 
of some off-road routes. These mapping tools do, however, 
when used in conjunction, provide an useful estimate of 
the total infrastructure footprint in this region. This is 
likely to be a conservative estimate due to our inability to 
map all infrastructure associated with oil and gas activity, 
or to account for indirect impacts from infrastructure 
presence such as those upon fauna such as the saiga. 

The direct footprint of the oil and gas sector is small 
compared with the area over which it is spread 
(~ 100,000 km 2 ), in fact constituting only approximately 
0.2% of the plateau by area, but this does not necessarily 
mean that impacts are negligible on the landscape scale. 
Measurement of the direct footprint does not account for 
the fact that infrastructure may well influence the ecology 
of a much wider area, for instance by changing vertebrate 
behavior up to 5 km away from infrastructure (Bem'tez- 
Lopez et al. 2010). The distribution of infrastructure is 
also important; as, beyond causing direct disturbance to 
wildlife, extensive linear infrastructure may physically 
fragment the landscape, which would particularly affect 
wide-ranging species (Rytwinski and Fahrig 2012) such as 
saigas. Quantifying these impacts is beyond the scope of 
this article, but these considerations are clearly important 
in designing biodiversity offsets, and in conservation 
planning for the wider landscape. 

Designing biodiversity offsets 

Biodiversity offsetting has been proposed as one mecha- 
nism within a raft of initiatives for reducing negative eco- 
logical impacts associated with infrastructure development 
on the Ustyurt (UNDP 2010). Offsets might enable the res- 
toration of highly degraded areas of Ustyurt habitat, com- 
mensurate with the land cleared for further infrastructure 
development (UNDP 2010; Bull et al. 2013b). A prerequi- 
site for implementing offsets is the evaluation of biodiver- 
sity losses, for which some metric based on both condition 
and area of habitat impacted is often used (Quetier and 
Lavorel 2011). Our estimates of total disturbance at land- 
scape level give a basic condition-area metric that could be 
used in partially quantifying the vegetation offset require- 
ments for the region. The estimate relates to development 
that predates any national offset policy, and sums habitat 
lost over several decades, so does not constitute actual off- 
set requirements. Furthermore, counterfactual trends in the 
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landscape, including the influence of other industries, 
human settlement and agriculture, and the influence of the 
Aral Sea crisis (Bull et al. 2013b), would need to be consid- 
ered in the design of any offset policy for the plateau. How- 
ever, our approach does outline a basic methodology for 
calculating vegetation loss in the region and contextualizes 
the magnitude of potential additional losses. 

More generally, our study shows that offset projects 
should consider linear infrastructure as well as hub devel- 
opment sites themselves, as linear infrastructure can consti- 
tute a large proportion of the total area impacted by 
developments. Finally, this investigation highlights the 
assumptions that must be made in calculating even one 
type of ecological loss for biodiversity offset schemes. In 
turn, this supports careful consideration of uncertainty in 
the development of offset policies, and the need to set up 
mechanisms to account for this uncertainty such as "offset 
multipliers" (Moilanen et al. 2009) and conservation bank- 
ing (Bekessy et al. 2010). In the development of offset pol- 
icy, it can often be assumed that the uncertainty dealt with 
by multipliers and biodiversity banking is that associated 
with the offset action (the biodiversity gain from the off- 
set), and not the quantification of the original impact (e.g. 
Defra 2011). But Moilanen et al. (2009) and Bekessy et al. 
(2010) suggest instead that multipliers include the uncer- 
tainty around the damage incurred, that is, the conserva- 
tion value of land that is developed, and that biodiversity 
banking should recognize uncertainty in impacts as well as 
actions. Considering uncertainty in the quantification of 
residual development impacts on biodiversity, as is carried 
out in our preliminary study here, is consequently crucial 
for the implementation of offset actions. 
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